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ABSTRACT 


dhk, 


The electric current helicity density x = {eijpbi ——) contains six terms, where hi are 

components of the magnetic field. Due to the observational limitations, only four of 
the above six terms can be inferred from solar photospheric vector magnetograms. 
By comparing the results for simulation we distinguished the statistical difference of 
above six terms for isotropic and anisotropic cases. We estimated the relative degree 
of anisotropy for three typical active regions and found that it is of order 0.8 which 
means the assumption of local isotropy for the observable current helicity density 
terms is generally not satisfied for solar active regions. Upon studies of the statistical 
properties of the anisotropy of magnetic field of solar active regions with latitudes and 
with evolution in the solar cycle, we conclude that the consistency of that assumption 
of local homogeneity and isotropy requires further analysis in the light of our findings. 

Key words: Sun: activity - Sun: magnetic field - Sun: Helicity 


1 INTRODUCTION 

Solar magnetic cycle is believed to be excited by solar dy¬ 
namo mechanism based on a joint action of differential rota¬ 
tion and mirror-asymmetric convection. Differential rotation 
known from helioseismology produces toroidal large-scale 
magnetic field from poloidal one while mirror asymmetric 
convection is responsible for transformation of the toroidal 
large-scale magnetic field into poloidal one in order to close 
the chain of self-excitation. Simple symmetry arguments 
show that a link between the toroidal and poloidal magnetic 
fields must be governed by a mirror-asymmetric quantity 
(generally a pseudo tensor) , however, there are sev¬ 
eral ways on how to implement this link in particulars (e.g. 
by notion of cyclonic motions, Parker (1955); or by mag¬ 
netic tubes arising and being twisted by Coriolis force, as in 
Babcock (1961), Leighton (1969) mechanism, what gives a 
variety of solar dynamo models. From the other hand, the 
degree of mirror asymmetry is believed to be moderate and 
the very degree of mirror asymmetry seems to be hardly de¬ 
termined from observational data. The point is that in order 
to quantify the mirror asymmetry of a dynamo one has to 
know 3D distribution of the mean-field characteristics and 
compute their spatial derivatives. 
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For example, a more straightforward quantity known 
as hydrodynamic (or kinetic) helicity density (v • curlv) 
which determines the excess of right-hand helixes against 
left-hand ones requires 3D distribution of velocity field v 
and its derivatives while conventional Doppler-effect gives 
line-of-sight velocities only. Note that we are interested in 
averaging quantities and (...) means corresponding aver¬ 
aging. While in mean-field dynamo theory the averaging is 
carried out over the ensembles of turbulent pulsations, prac¬ 
tically given the observational data in vector magnetograms 
of solar active regions, we average by area in the available 
field-of-view of something similar to that. 


A practically accessible way to quantify mirror asymme¬ 
try observationally was firstly suggested by Seehafer (1990). 
He pointed out that vector magnetographic observations 
provide three magnetic field components on a surface z = 
const at solar atmosphere (local coordinates x and y are 
parallel to the solar surface, and for a limited field of view 
we can ignore the effect of curvature). He used the force-free 
field parameter a as a proxy of electric current helicity. The 
available vector magnetic field is not sufficient to calculate 
the entire current helicity density x = (H • curlB) and to 
quantify relative number of right handed twisted magnetic 
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tubes against left-handed ones, however, one can calculate 
a mirror-asymmetric quantity 

= (B^(curlB)^) = {Bz {VxBy - VyB^)) , (1) 

which makes one of the three believed to be similar con¬ 
tributions to the entire quantity y, i.e., Xcc and Xy (fhey 
can be obtained from Eq. (1) by circular replacement of in¬ 
dices). Note that Xz do not contain derivative in direction z. 
If magnetic field B is locally statistically isotropic, all three 
contribution in x are statistically equal and the natural con¬ 
sequence is 

X = ^Xz- ( 2 ) 

The force-free field parameter a and the cur¬ 
rent helicity parameter x become available for observa¬ 
tional determination in solar active regions (see, e.g., 
Pevtsov 1994; Abramenko et al. 1996; Bao & Zhang 1998; 
Hagino and Sakurai 2004). The above statistical studies 
show that the sign of a and x is same, which is predom¬ 
inantly negative in the northern hemisphere and positive in 
the southern hemisphere. Results of monitoring x^ in two 
last solar cycles and butterfly diagrams for solar cyclic vari¬ 
ation of this quantity are presented in Zhang et al. (2010). 
The quantity followed the helicity polarity rule as well as 
pronounced areas on the butterfly diagrams where the po¬ 
larity rule is inverted (Zhang et al. 2010). The result looks 
instructive for solar dynamo modeling (Zhang et al. 2012). 
Dynamo interpretation of the observational data is usually 
based on assumption of the local statistical isotropy, and so 
Xz is considered as an observational tracer for x- 

Conventional theory of turbulence is originated from 
Kolmogorov (1941), while Iroshnikov (1963) and Kraichnan 
(1965) ideas presumes that the velocity and magnetic fields 
become statistically homogeneous in sufficiently small scales 
(though Goldreich and Sridhar (1997) have stressed the role 
of magnetic field anisotropy in MHD turbulence). It is how¬ 
ever not clear what scales in the solar photospheres can be 
considered as homogeneous enough for this assumption and 
whether typical active regions may fall under this consider¬ 
ation. A perspective to verify to what extent the solar mag¬ 
netic field at the scale of active regions can be considered 
as statistically isotropic has not been yet considered obser- 
vationally at least in the approach suggested by Seehafer 
(1990): the point is that calculation of other contributions 
to the current helicity to be compared with Xz would require 
derivatives of magnetic field components in z direction while 
in fact we have three components of the magnetic field at 
surface z = const only. 

The aim of this paper is to show that one can refur¬ 
nish the approach of Seehafer (1990) in a way to verify the 
hypothesis of local isotropy in the scales of active regions. 
Instead of to present x as a sum of three contributions x^? 
Xx and Xy which have to be equal in a locally isotropic case, 
we present this quantity as a sum of six quantities to be 
equal in the isotropic case. Four out of the six quantities 
do not contain derivatives in z direction and are potentially 
accessible for observations. 

We demonstrate that the available bulk of data for the 
magnetic field vector in solar active regions enable us to 
obtain statistically robust estimates for the four quantities. 
Local statistical isotropy implies two pairs of identities for 


the above four quantities. We show that one pair of the 
identities holds while the other fails being confronted with 
the observational data. We may believe that the magnetic 
field occurs to be substantially statistically anisotropic in 
scales of active regions. As a result, the quantity x^ which 
has been traced for last two solar cycles has to be considered 
as a specific mirror asymmetric tracer of anisotropic solar 
MHD rather than purely x/3. We discuss the importance of 
this conclusion for solar dynamo models. 


2 OBSERVATIONAL DATA 

We used 6629 vector magnetograms observed by Solar Mag¬ 
netic Field Telescope (SMFT) at Huairou Solar Observing 
Station from 1988 to 2005. This data sample was used earlier 
by Gao et al. (2008) and Zhang et al. (2010). The SMFT is 
equipped with a birefringent filter for wavelength selection 
and KD*P crystals to modulate polarization signals. The 
Fe 1 5324.19 A line is used. A vector magnetogram is built 
using four narrow-band (0.125A) filtergrams of Stokes I, Q, 
U and V parameters. The center wavelength of the filter 
can be shifted and is normally at -0.075 A for the measure¬ 
ments of longitudinal magnetic field and at the line center 
for the transversal magnetic fields (Ai & Hu 1986). The 180° 
ambiguity in the azimuth angle (0) was resolved following 
Wang, Xu, and Zhang (1994) by comparison with a poten¬ 
tial field. We used the method given by Gao et al. (2008) 
to correct the Faraday rotation to the azimuthal angles for 
vector magnetic field. 

We have also established the levels of noise for the lon¬ 
gitudinal and transverse components of the magnetic field 
as 20 G and 100 G, in accord with previous works on anal¬ 
ysis of Huairou SMFT data (e.g., Abramenko et al. 1996). 
That values have been used to estimate the impact of noise 
in statistical studies as it has been performed earlier by e.g. 
Bao & Zhang (1998) or Zhang et al. (2010). 


3 THE METHOD 

By definition the current helicity is a scalar product of the 
magnetic field pseudo-vector and its curl (proportional to 
the electric current) vector. Generally speaking this con¬ 
struction can be obtained as a trace of a more generic 
three-index tensor quantity as it has been considered by 
Zhang et al. (2012), see their Appendix 1. For our analy¬ 
sis, let us denote local observable quantities by low cases, 
i.e. the local magnetic field b which is observable at the so¬ 
lar surface z = 0 as well as the local electric current helicity 
density /ic = b • curlb. We denote (x, y) local Gartesian co¬ 
ordinates on the image plane. The direction to the observer 
z on the plane is fixed z = 0. We reserve above mentioned 
notations B and x foi* fhe magnetic field and current he¬ 
licity in homogeneous and isotropic model. Then according 
to definition of curl this quantity naturally comprises of six 

parts: x = = hi + /i 2 + /13 + where 
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Figure 1. The vector magnetograms for three examples: NOAA 
8898 (S13.0W7.0) observed by SMFT on 01:34 UT, March 8, 2000 
(top); NOAA 6659 (N28.6W4.5) observed by SMFT on 05:29 UT, 
June 9, 1991 (middle); NOAA 11158 (S19.0E11.0) observed by 
HMI on 23:59 UT, Feb 14, 2011 (bottom). The arrows represent 
the direction of transverse field and the contours represent the 
longitudinal magnetic fields of ±100, 400,1600, 3200G. 


(%) ’ '‘I - *• ( 


We denote integral quantities over the available magne- 
togram field-of-view by capital cases Hi — J hidxdy for 
z = lto6, so the overall average current helicity reads 


He = iFi±iF2±i^3±i^4±i^5±i^6 (4) 

+ J bx dxdy+ J bx 
+ J by (^^^ dxdy+ J by (^“^^ dxdy. 

Using the integral by parts formula, we can obtain the 
following relation between helicity parts, for example. 



bz hydl 


(5) 


where the latter integral is taken over contour F at the 
boundary of our field of view (boundary integral). If we as¬ 
sume the magnetic field at the boundary of the active regions 
is weak enough, that means the boundary integral is very 
small, then the averages of the corresponding parts helicity 
are approximately equal in pairs Hi ~ He and H 2 — H 3 . 
We illustrate validity of this consideration below. 

For the purpose of illustration of our data we selected 
two active regions observed by SMFT in 1999 and 2000, and 
one to compare with other instruments, we used the data for 
an active region observed by Helioseismic and Magnetic Im¬ 
ager onboard the Solar Dynamics Observatory (HMI/SDO) 
in 2011. Figure 1 shows the vector magnetograms for three 
active regions observed by SMFT and HMI/SDO. We can 
see that the magnetic field is nearly potential for NOAA 
8898 and the magnetic field is strong helical for NOAA 6659. 
The spatial resolution of HMI is higher than SMFT. We 
calculated the integrals Hi and He and the boundary inte¬ 
grals in formula (5). The results are listed in Table 1. The 
boundary integral is about 2.87%, 0.31% and 3.48% of the 
mean value of and He for active region NOAA 8898, 6659 
and 11158 respectively. In order to estimate tolerable error 
in computation of the boundary integral we use the typical 
noise levels: 20 G in Bz and 100 G in For example, active 
region NOAA 6659, the error is about 2.04% of the mean 
value of and He- The corresponding boundary integral 
is well below the error level. These observational examples 
show that the above assumption of validity of integration 
by parts works pretty well when the magnetic field at the 
boundaries of the field of view is weak, so the formula (5) 
of integration by parts results in these equalities. Therefore, 
we established on the basis of both theoretical consideration 
and observational illustration that the four parts of observ¬ 
able current helicity are equal by pairs to their counterparts: 
Hi ^ He and H 2 ^ H 3 . 

It is trivial to see that the local values of ah helicity 
parts are generally unequal. Now let us consider whether 
the integral identities implied to mean helicity parts hold 
for the observational data. 

First of ah, ah the parts of the total current helicity 
mentioned in the Introduction can be expressed in notations 
analogous to the ones in formulae (5) and their observational 
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Table 1. Summary of illustration computations of the helicity and boundary integrals in formula (5). The unit of Hi, Hq and boundary 
integral is 


Instr. 

NOAA 

date 

Time 

coordinate 

Hi 

He 

% difference 

bound.int. 

err.est. 

SMFT 

8898 

2000.03.08 

01:37 

S13.0W7.0 

-0.0493 

-0.0495 

2.87% 

0.0014 

0.0135 

SMFT 

6659 

1991.06.09 

05:29 

N28.6E4.5 

-1.4978 

-1.4983 

0.31% 

-0.0046 

0.0306 

HMI 

11158 

2011.02.14 

23:47 

S20W17 

0.0793 

0.0784 

3.48% 

-0.0027 

0.0123 


- Xl - /b2 - X3 X4 - X5 - T6 




-0.5 0.0 0.5 1.0 1.5 


Figure 2. PDF for the six helicity parts: non-helical isotropic 
case (top), helical isotropic case (bottom). Note that only four 
out of the six are of observational interest. The x-axis is the rela¬ 
tive value of helicity parts and the y-axis is the probability density 
function. All parts statistically coincide with each other though 
for the helical case there is a pronounceable bias visible as asym¬ 
metry over the center of distribution. 

counterparts 

X-s — Xi + X2 ^^ Hi + H2 = Hcz , 

Xcc = X3 + X4 < - > H3 + H4 = Hex , (6) 

Xy — X5 + X6 ^ H^ He = Hey , 

and notice that only the former part can be fully computed 
from observations as it does not contain derivatives with 
respect to z. Only one term in each of the two latter parts 

can be computed and the two terms H 4 and He are not 

available from observations on the image plane. We either 
cannot use formula (5) for them as the relevant derivative is 
with respect to z but integration is carried out over x and 

y- 


Let us assume local statistical isotropy of turbulence 
and take for instance two parts of helicity equal, say 

X^ = Xx 

Then we immediately have that as iLi + H 2 = H 3 + H 4 , 
and due to formula (5) Hi = He and H 2 — H 3 , then the 
three parts are equal 

Hi = H4 = He. (7) 

Therefore, an additional assumption on the equality of the 
other helicity parts 

Xx = Xy 

would automatically lead to equality of the other three 
parts 

H2 = H3 = He. ( 8 ) 

The above consideration means that for veriheation of 
the assumption of local isotropy unobservable parts of he- 
licity Ha and He need to be evaluated in order to check 
equations (7-8). 

Before doing further analysis we are going to see what 
relationship between the parts of helicity we can expect from 
theoretical consideration. For that purpose model simulation 
of the helicity parts is performed in the next section. 


4 SIMULATIONS 

To prescribe a quasi-random magnetic held B with vanishing 
mean value in a periodic box, we use a Fourier expansion 
in modes with randomly chosen directions of wave vectors 
k but with amplitudes adjusted to reproduce any desired 
energy spectrum: 

/ B(k)e‘'' "d^k, (9) 

where B is the Fourier transform of B. The corresponding 
magnetic energy spectrum is given by 

Mik)= [ |B(k')|"dV, (10) 

J\\<.'\=k 

where the integral is taken over the spherical surface of ra¬ 
dius k in /c-space. In the isotropic case, M(k) = 47r/c^|B(A:)|^. 
In order to ensure periodicity within a computational box 
of size L, as required for the discrete Fourier transforma¬ 
tion, the components of the wave vectors are restricted to 
be integer multiples of 27r/L. 
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Figure 3. PDF for the six helicity parts same as in Figure 2: non-helical anisotropic case for anisotropy in a:-direction (top left) and 
^-direction (top right), helical anisotropic case for anisotropy in x-direction (bottom left) and z-direction (bottom right). 


Table 2. Summary of the model simulation cases partially presented on Figure 3. Note that observationally available parts of helicity 
are approximately equal in the case with anisotropy with respect to 2 ; only. They are equal in all helical isotropic cases while in other 
anisotropic helical cases they are not equal. 


case (mean ± dispersion) 

XI 

X2 

X3 

X4 

X5 

X6 

Isotopic Non-helical 

0.A2.71 

0.A2.74 

0.A2.62 

0.A2.61 

0.A2.67 

0.A2.67 

Isotropic Helical 

1.56A3.12 

1.53A3.12 

1.53A3.11 

1.55A3.12 

1.56A3.12 

1.56A3.11 

Anisotropic-x Non-helical 
Anisotropic-y Non-helical 
Anisotropic- 2 : Non-helical 

0.A0.19 

0.A1.81 

0.A2.2 

0.A1.93 

0.A0.2 

0.A2.13 

0.A2.28 

0.A0.2 

0.A1.9 

0.A2.13 

0.A1.88 

0.A0.2 

0.A1.79 

0.A2.08 

0.A0.2 

0.A0.19 

0.A2.14 

0.A1.82 

Anisotropic-a: Helical 
Anisotropic-y Helical 
Anisotropic- 2 ; Helical 

0.05A0.2 

1.29A2.35 

1.26A2.44 

1.29A2.24 

0.07A0.22 

1.26A2.47 

1.29A2.41 

0.07A0.22 

1.26A2.19 

1.36A2.57 

1.38A2.34 

0.06A0.2 

1.23A2.04 

1.25A2.44 

0.06A0.2 

0.05A0.2 

1.29A2.54 

1.26A2.18 


A solenoidal vector field B, i.e., that having k-B(k) = 0, 
is specified by 




Random choice a complex vector X implies zero net cur¬ 
rent helicity of B. We consider a magnetic energy spectrum 


represented by two power-law ranges, 


M{k) = Mo 


{ 




for k < ko, 
for k > ko, 


( 11 ) 


with So > 0, Si > 0 and Mq = 1, where /co = 6 is the 
energy-range wave-number. We use si = 5/3 as in Kol¬ 
mogorov’s spectrum (Kolmogorov 1941) and so = 2 as in 
Christensson et al. (2001). 
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Figure 4. Distribution of x^:, xi cind X 2 for three cases: non-helical magnetic field (top), helical magnetic field (middle), helical magnetic 
field plus potential magnetic field (bottom). 


The helical B can be obtained with choice X as 


X = Yzbi|Y| 


k X Y 

|k X Y|’ 


( 12 ) 


where the sign defines the sign of current helicity 
and Y(k) = Y(—k) is a random real vector. Condition 
(12) implies B(k) = B*(—k). 


The probability density functions (PDFs) calculated 
from 2D distributions [(x,y)-plane of 3D simulated cube] 
are shown in Fig. 2. First of all we note that isotropy means 
similarity of distributions of xi? X 2 , Xs cind xe- Secondly, 
non-zero helicity leads to asymmetry of PDFs. Furthermore, 
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Figure 5. The examples of PDFs for the four observable helicity parts: a) column is for NOAA 8898; b) column is for NOAA 6659; c) 
column is for NOAA 11158. The upper (bottom) panels are the PDFs for helicity parts hi and Hq (h 2 and hs). The mean values and 
the standard deviations of helicity parts are given on the top right corner of each panel. 


we can produce corresponding statistically anisotropic fields 
for non-helical and helical cases. 

Anisotropic case is simulated by the additional factor 
in M{k). One can take 


M(k) = k~^Mo 


{ 


{k/koyo 


for k < ko, 
for k > ko, 


(13) 


The corresponding PDFs are shown in Figure 3. From Fig¬ 
ure 3 one can see that anisotropy in different direction af¬ 
fects the PDFs of different parts of helicity in a different way. 
Anisotropy in x-direction on the left two panels leads to the 
parts of helicity xi and xe containing derivatives in this di¬ 
rection to be distributed with much lower dispersion. One 
can see the same for anisotropy in ^-direction in the right 
panels for X 4 and xs- Each of the two quantities in these 
pairs are distributed statistically similar. For the remaining 
parts of helicity, in the helical case (bottom panels), we can 
see that the left and right tails of the PDFs have different 
spreads and in particular the left tails are more inclined. 
Furthermore, we see that these four remaining parts of he¬ 
licity group in two statistically similar pairs, with respect to 
the component of the magnetic field which enters into each 
of the parts, namely Bx for xs cind X 4 in case of anisotropy 
in X— direction and Bz for xi and X 2 in case of anisotropy 
in ^—direction. 

The results for the cases shown in Figure 3 are summa¬ 
rized in Table 2. One can see that the mean values of the six 
helicity parts are approximately equal in the isotropic heli¬ 
cal case. In the anisotropic helical cases, in accord with the 
integration by parts formula (5), the mean values of helicity 
parts are equal in pairs: (xi) ~ (xe) and (x 2 ) ~ (xs)- The 
other two parts (x 4 ) and (xs) are not subject of the formula 
(5) as they contain derivatives with respect to z while av¬ 
eraging involves differentiation over the image plane (x,y). 


and they may generally be different, as it is in the anisotropic 
helical case. 

Now let us consider how the sign of the parts xi? X 2 , 
X 3 and X6 can locally represent the sign of the total helicity. 
One effect can be the contribution from the potential mag¬ 
netic field. We simulate distributions of x^? Xi and X 2 for 
three cases: non-helical magnetic field, helical magnetic field 
and helical magnetic field plus potential magnetic field (see 
Figure 4 from top to bottom). As expected for the purely 
non-helical magnetic field, the all three maps possess the 
same kind of patterns with alternating sign. For the heli¬ 
cal case we also have the same kind of patterns but with 
dominating sign of helicity. Additional contribution from the 
potential magnetic field does not change much the total he¬ 
licity x^ but its parts xi and X 2 have strongly alternating 
pattern unlike for the case of non-helical magnetic field. We 
can see that assigning these two parts the opposite signs in 
the alternating pattern may merely cancel each other. 

Therefore, we have seen that the four parts of helicity 
that have observational interest are close in pairs, but be¬ 
tween the pairs there may be a significant difference in their 
distributions and the integral values. We have also noted 
now the specific contribution of one and the other parts may 
partially cancel each other in the overall helicity. 


5 OBSERVATIONAL RESULTS 

Now we compare theoretical predictions with results of ob¬ 
servational data analysis. Figures 5 show PDFs for the four 
observationally available parts of helicity hi, ^ 2 , ha, he com¬ 
puted for all pixels in the magnetograms of the three active 
regions shown in Figures 1. One can see that the mean val¬ 
ues of distribution for pair (hi, he), and pair (h 2 , hs) are 
very close to each other with accuracy of a few per cent. 
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Figure 6. Distribution of hi and /12 for the three examples: Nearly potential magnetic field in NOAA 8898 (top), helical magnetic 
field in NOAA 6659 (middle), relatively high special resolution magnetic field in NOAA 11158. The unit of helicity is 10“^G^/m. 


but the difference between the mean values of hi and /12 {h^ 
and ho) is large although the distribution pattern is similar. 
This observational result is similar with the anisotropic case 
simulation in Figure 3 although some difference is found. 
This discrepancy can be attributed to the complexity of the 
observation. We also show the distribution map for hi 
and /i 2 in Figure 6. We can see that there is a significant 
difference in the distribution of hi and /i 2 , but their sum is 
nearly equal to hz- 

We have analyzed similar distributions for various vec¬ 
tor magnetograms obtained at Huairou Solar Observing Sta¬ 
tion (HSOS) available for date in cycles 22 and 23 as well as 
a recent vector magnetograms observed by HMI/SDO in or¬ 
der to learn that the properties above for this PDF is quite 
generic for the magnetograms under consideration. 

Now from the study of distribution of the parts of helic¬ 
ity by pixels over a given magnetogram of an active region 
we move towards study the distribution of the mean val¬ 
ues of these quantities for 6629 magnetograms observed by 
HSOS from 1988 to 2005. 

Figure 7 shows scatter plots for helicity parts for our 
sample. One can see that correlation between Hi and Hq 


{H 2 and is very high while the one between Hi and H 2 
{H 3 and He) is low. This means that average with signed 


flux bz weighted by 


dbx 

dy 


is systematically different from 


the one weighted by —-— in a magnetogram. These con- 
ax 

Arm that the properties which we see above for the PDFs 
of helicity parts in a magnetogram: robustness of numeri¬ 
cal scheme for the use of integration by parts, as well as 
absence of isotropy. Comparing with our theoretical simu¬ 
lation in the section above, this statistical observational re¬ 
sult corresponds to anisotropic helical case while anisotropy 
is either by x or y, or both. This is not the case of purely 
anisotropic by only vertical ^-direction (convective stratified 
turbulence). We also studied the sign agreement between he- 
licity parts. We found that 95.5% (95.9%) of Hi and Hq {H 2 
and H 3 ) among 6629 vector magnetograms agree in their 
signs, while for Hi and H 2 {H 3 and Hq) the percent of agree¬ 
ment is much lower which is 38.0% (34.7%). We checked the 
vector magnetograms whose Hi and Hq {H 2 and H3) have 
different sign and found that the magnetic fields are not 
weak at the boundary for some of those magnetograms. On 
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Figure 7. Scatter plots key helicity parts over the dataset 1988- 
2005. each point represents average value of helicity parts for one 
magnetogram. The unit of helicity is 10“^G^/m. 


the other hand, the values of helicity parts inferred from 
those magnetograms are often smaller by one or two orders 
of magnitude than for the typical magnetograms. Those val¬ 
ues are too small to be estimated accurately. 

Let us study how the local anisotropy of current helicity 
varies with the solar cycle. From Figure 8, we can see that 
the latitudinal structure and variation with the solar cycle 
of {Hi) and {Hq) ( {H 2 ) and {H 3 )) are identical except a 
small part which is different. The reasons for this difference 
are analyzed in the above paragraph. While the latitudinal 
structure and variation with the solar cycle of {Hi) and {H 2 ) 
( {H 3 ) and {He)) are much different which conhrms again 
the absence of isotropy for the observable current helicity in 
solar active regions. We can also see that in some periods 
and latitudinal intervals the helicity parts satisfy the hemi¬ 
spheric sign rule (negative/positive sign in the North/South 
hemispheres) well, e.g., 1992-1994 and 2001-2003 for {Hi) 
in southern hemisphere. The fractions of magnetograms fol¬ 
lowing the hemispheric sign rule are 46.2% for {Hi), 61.9% 
for {H 2 ), 60.2% for {H 3 ) and 47.2% for {He) in northern 
hemisphere respectively, while in southern hemisphere they 
are 57.9% for {Hi), 55.4% for {H 2 ), 54.4% for {H 3 ) and 
57.1% for {He) respectively. 

The sum {Hcz) = (i7i) + {H 2 ) usually used in the past 
study in Figure 9 shows some anti-symmetry with respect 
to the solar equator and the hemispheric sign rule is well 
pronounced, this hgure is almost the same as Figure 2 in 
Zhang et al. (2010). The use of the sum of the other two he- 
licity parts He and H 3 , respectively, would produce visually 
identical result, and not shown here. The fractions of mag¬ 
netograms following the hemispheric sign rule for {Hcz) are 
58.9% in northern hemisphere and 61.7% in southern hemi¬ 
sphere respectively. We can see some regular inversions of 
hemispheric sign rule within isolate ranges of latitudes near 


the beginning and the end of the solar cycle for {Hcz)- We 
note that during these phases the signs of pair of helicity 
parts {{Hi) and {H 2 )) are often the same and violate the 
hemispheric sign rule, which makes their joint contribution 
to the reversal of hemispheric sign rule. 


6 DISCUSSION 

In order to establish the degree of anisotropy quantitatively, 
we take 3 active regions which used for calculating boundary 
integral in section 3 for examples. We estimate the relative 
degree of anisotropy, for example, as a relative anisotropy 
imbalance between the two parts of helicity using some norm 
as following: 


qi2 


l|4i - h2\\ 
l|/il|l + l|42||’ 


where the norm ||/i|| can be computed as 


(14) 


where w is weight factor as follow: when we integrate in the 
magnetogram using ah pixels, then re = 1. Alternatively, we 
can use only pixels where the signal is greater than cut-off 
noise levels {\Bz\ > 20G and Bt > lOOC) as we usually do 
so in helicity statistical studies (re = 1 if greater, or re = 0 
if less than noise level). The value of qi 2 is between 0 to 
1. The results are listed in Table 3. qi 2 is calculated using 
ah pixels in the magnetogram, qi 2 is only using the pixels 
where the signal is greater than cut-off noise levels. From 
this table, the order of qi 2 is from 0.74 to 0.83 and it is not 
affected much by the noise, which is close to 1, it says that 
the observed quantities are extremely anisotropic. 

We may further speculate about the sources of 
anisotropy of current helicity. One effect which sounds trivial 
in the solar photosphere would be vertical stratihcation of 
convection. We, however, are inclined to discuss this issue as 
the parts of helicity involving vertical derivatives of the mag¬ 
netic held are not observationahy available. The difference 
in properties and behaviour of i7i and H 2 is likely a mani¬ 
festation of the effect of rotation as one contains derivatives 
with respect to azimuthal (cohinear to rotation), and the 
other meridional direction (perpendicular to rotation). The 
other source of anisotropy could be the effect of large-scale 
magnetic held. Two components (azimuthal and radial) of 
the solar cyclic magnetic held are anti-symmetric over the 
equator and the other meridional component is symmetric 
over the solar equator. 

The complete use of the assumption of local statistical 
isotropy means also the equality between these two groups 
of the three parts of helicity in equation 7 and 8. Two of 
them, Ha and He, however, can not be calculated from the 
observational vector magnetograph data at the solar photo¬ 
sphere. Therefore, we can not check the condition of local 
isotropy of solar turbulence in active region by observation¬ 
ahy available data completely. On the Other hand, there 
are several factors which affect the precision of the calcu¬ 
lated helicity parts, such as that the accuracy of transversal 
held is lower than the one of longitudinal held, the magneto¬ 
optical effects, calibration of magnetic helds etc. We cannot 
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1988 1994 2000 2006 

Time (year) 

• •••ooOQ 

- 0.43 - 0.43 — 0 . 28 - 0 . 28 — 0 . 12 - 0 . 12 - 0 0 - 0.26 0 . 26 - 0.42 0 . 42 - 0.57 > 0^7 



1988 1994 2000 

Time (year) 

• • • • o O O 

<- 1.04 - 1 . 04 — 0 . 70 - 0 . 70 — 0 . 35 - 0 . 35 - 0 0 - 0.15 0 . 15 - 0.36 0 . 36 - 0.57 


2006 


Figure 8. Evolution of four helicity parts with the solar cycle overlaid with sunspot density (color). The vertical axis gives the latitude 
and the horizontal axis gives the time in years. The circle size gives the magnitude of helicity as averaged over two-year running windows 
over latitudinal bins of 7° wide. The unit of helicity is 


Table 3. Summary of the relative degree of anisotropy. The unit of Hi and H 2 is 10^'^G^m. 


Instr. 

NOAA 

date 

Time 

coordinate 

Hi 

H 2 

qi 2 

9f2 

SMFT 

8898 

2000.03.08 

01:37 

S13.0W7.0 

-0.0493 

0.0310 

0.8348 

0.8325 

SMFT 

6659 

1991.06.09 

05:29 

N28.6E4.5 

-1.4978 

-2.5897 

0.7823 

0.7794 

HMI 

11158 

2011.02.14 

23:47 

S20W17 

0.0793 

0.7739 

0.7472 

0.7472 



Figure 9. Evolution of the observationally available helicity de¬ 
fined as a sum of two parts + {H 2 ) with the solar cycle. The 
unit of helicity is I'm. 


completely estimate the influence of these factors on the he¬ 
licity parts although we used some data reduction method 
to reduce it. This topic requires further investigation. 


7 CONCLUSION 

We have studied the distribution and properties of the parts 
of current helicity from observation and simulation. The 
main conclusions are as following: 

(1) The simulation results show that the means of the 
six helicity parts statistically coincide with each other for 
the isotropic case, while for the anisotropic helical case it is 
different. 

(2) The distribution of the observed helicity parts is 
similar with the helical anisotropic case in simulations. This 
shows that the distribution of the observed helicity parts is 
anisotropic over active region scales. Assumptions of local 
homogeneity and isotropy in computation of observational 
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proxies of helicity require further analysis in the light of our 
findings. 

(3) The four observable helicity parts are equal in pairs 
owing to the magnetic field is weak at the boundaries of 
magnetogram, but there is a large difference between dif¬ 
ferent pairs which may be caused by the anisotropy of the 
current helicity density in active region. Both the pairs of 
the helicity parts follow the hemispheric sign rule with cer¬ 
tain exceptions in periods and hemispheres but the sums of 
the two parts follow the hemispheric sign rule more robust 
and uniform than the individual parts of helicity alone. 

More theoretical modeling of anisotropy in solar-like 
turbulence is required to understand these statistical results. 
Our simple simulations have shown that the anisotropy may 
be present in several directions and not only in the direction 
of vertical stratification. This is confirmed by several exam¬ 
ples of observational data as well as their statistical anal¬ 
ysis. Other sources of anisotropy such as inhomogeneity or 
rotation, and furthermore presence of the large-scale mean 
magnetic field which alternates in sign with every 11-year 
sunspot cycle may be possible explanations of these results. 
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